Soliton acceleration by dispersive radiation: a contribution to rogue waves? 
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Rogue waves are solitary waves with extreme amplitudes, which appear to be a ubiquitous phe- 
nomenon in nonlinear wave propagation, with the requirement for a nonlinearity being their only 
unifying characteristics. While many mechanisms have been demonstrated to explain the appear- 
ance of rogue waves in a specific system, there is no known generic mechanism or general set of 
criteria shown to rule their appearance. Presupposing only the existence of a nonlinear Schrodinger- 
type equation together with a concave dispersion profile around a zero dispersion wavelength we 
demonstrate that solitons may experience acceleration and strong reshaping due to the interaction 
with continuum radiation, giving rise to extreme-value phenomena. The mechanism is independent 
of the optical Raman effect. A strong increase of the peak power is accompanied by a mild increase 
of the pulse energy and carrier frequency, whereas the photon number of the soliton remains prac- 
tically constant. This reshaping mechanism is particularly robust and may explain the appearance 
of rogue waves in a large class of systems. 

PACS numbers: 42.65.-k, 42.81. Dp, 47.35.Fg, 04.70.-s 



The appearance of waves with extreme amplitude has 
been investigated in a large class of physical systems 
[1-8 . Their appearance is most drastically illustrated 
for the case of ocean waves P4TT], with waves exceed- 
ing the average wave crest by a factor three or more 
and causing serious damage to ocean-going ships. Re- 
cently, similar phenomena have also been reported in op- 
tics in the soliton-supporting red tail of supercontinua 
(SC) in fibers [8 . Substantial progress has been made 
in the understanding of the mechanisms behind optical 
rogue waves p!2UT5] . Currently, most explanations fol- 
low one of two alternatives. One involves soliton fission 
and selective Raman shifting of the largest solitons to the 
long- wavelength side of the spectrum [131 E] • The other 
builds on the dynamics of Akhmediev breathers [K and 
inelastic collisions between solitons or breathers. While 
these explanations concentrate either on soliton fission or 
fusion processes, we demonstrate in the following that a 
similar mechanism may exist between solitons and con- 
tinuum radiation in the normal dispersion range. We will 
refer to the latter as dispersive wave (DW) . Suitable con- 
ditions provided, this DW can strongly modify a soliton 
through cross phase modulation (XPM). Strong reshap- 
ing, in particular temporal compression, of the soliton is 
accompanied by a mild increase of its energy and carrier 
frequency, while the photon number of the soliton re- 



mains practically constant. Decrease of the center wave- 
length gives rise to acceleration of the soliton by virtue 
of dispersion. The soliton peak power grows surprisingly 
and may more than double in only a few centimeters of 
propagation. 

In the following, we consider the rogue wave formation 
in the SC generation in a single-mode photonic crystal 
fiber with one zero dispersion wavelength (ZDW), which 
is similar to the fiber used in [SI [1^. In the following, 
we restrict our analysis to the minimum set of optical 
effects necessary for the formation of rogue waves. These 
effects include dispersion and the Kerr nonlinearity, yet 
exclude Raman scattering, cf. [16 . In the fiber geometry, 
the optical field is characterized by a single real valued 
component E{z^t) whereas dependencies perpendicular 
to the propagation coordinate z are integrated out. We 
choose a suitable time period T, introduce frequencies 
u e 27rZ/T, denote spectral field components by E^{z), 
and following [17 define a complex valued E{z^t) such 
that Suj{z) = E^{z) - idzE^{z)/\j3{uj)\. Note that E = 
Re[£^], the propagation equation for S{z^t) reads 

o 2 (3) 

id,£^ + \P{uj)\£^ + ^^^^Q£\^£)^ = 0. (1) 

Parameters c, x^^\ the speed of light, the 

third-order nonlinear susceptibility, and the propagation 
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FIG. 1: Temporal evolution of a higher- order soliton injected 
close to the ZDW into the fiber along z for a typical SC gener- 
ation process. Note that the calculation does not involve the 
Raman effect, (a) Final state with a rogue wave exceeding the 
average wave crest by more than a factor of three, (b) Prop- 
agation dynamics of the solitons and non-solitonic radiation, 
generated by the fission process. 



constant respectively. Equation ([T]) is subject to the con- 
servation laws 
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where n{uj) is refractive index and 1 1^2 are finite and pro- 
portional to the time-averaged photon flux and power re- 
spectively [iT. For unidirectional propagation E{z,t) is 
identical to analytic signal E{z^t) = 2 ^^^q ^a;(z)e~*^^ 
and moreover only the positive- frequency part of \E\'^E 
contributes to Eq. ([T]). The fiber propagation constant 
may be obtained by numerical integration of the group 
delay Pi{oj) = P'{oj) and then approximated following 
[18 . If the slow envelope description with respect to a 
carrier frequency ujq applies, Eq. ([T]) reduces to the stan- 
dard nonlinear Schrodinger equation [19] with the nonlin- 
earity parameter 7 = (3co'oX^^^)/[4eoc^^^(co'o)^eff]5 where 
^eff is the effective fiber area. Beyond the standard treat- 
ment with an envelope approximation, our approach cor- 
rectly models nonlinear processes between spectrally dis- 
parate waves, i.e., four- wave mixing processes and XPM 
between solitons and DWs, and between individual soli- 
tons. These nonlinear processes have been previously 
found important for explaining rogue waves in gravity 
matter waves [10]. 

We launch a hyperbolic secant pulse (center wave- 
length 897 nm, full width at half maximum FWHM = 265 
fs), corresponding to a higher-order soliton with soliton 
number A/" ^ 28 in the anomalous dispersion regime of 
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FIG. 2: (a) Visualization of the soliton propagating from A 
to B in Fig. [l] Temporal delays are shown relative to the 
unperturbed propagation of the soliton z — 1.5-2 cm. For 
clarity, the width of the soliton has been stretched by a factor 
5. Color coding visualizes Ao(^), which changes from 1060 
to 985 nm (red and blue, respectively), (b) Development of 
soliton parameters \q{z) (thick red line), pulse energy oc I2 
(solid black line), and peak power Po{z) (dashed line). Energy 
content of the dispersive wave within zbl.5T interval around 
t* [z) is shown as a dotted line. 



the fiber close to the ZDW = 842 nm. For a nonlinear 
fiber with 7 = 0.1W~^m~^ this corresponds to a peak 
power of 19 kW, so that soliton fission is favored [20] . 
These conditions ensure the formation of a SC with the 
increase of the initial spectral width by one to two or- 
ders of magnitude [21 . Figure [l] shows the typical SC 
evolution in the temporal domain. With the rather mod- 
erate powers in this example, the effect of the modu- 
lational instability can only be observed in the initial 
^ 0.5 cm propagation length before pulses reach the sub- 
100 fs regime [22 . The fundamental solitons produced 
in the fission process exhibit durations between 10 and 
20 fs with different peak powers, appearing as pronounced 
lines which clearly stand out from the background. The 
fission process also generates DWs in the normal disper- 
sion regime ^J. The further away from the ZDW soli- 
tons are being generated, the slower they will propagate 
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[22] , accumulating delay (Fig. [T]). As we deliberately ex- 
cluded Raman scattering in our analysis, we might expect 
that the group velocity of solitons is constant except in 
places where isolated soliton-soliton scattering processes 
occur. Indeed, inspection of Fig. [l] reveals several such 
characteristic crossings of soliton trajectories in the t-z 
plane. However, it also reveals that the trajectory of 
the strongest soliton does not appear to be ruled by rare 
isolated scattering events (AB, Fig. [l]). The parabolic 
trajectory of this soliton is witness of its constant accel- 
eration. 

To elucidate the physical mechanisms behind this pe- 
culiar acceleration, we numerically isolated the soliton, 
separated it from accompanying continuum radiation, 
and fitted the model function f{t) = Pq sech^[(t — t^)/to] 
to its intensity envelope (Fig. [2) FWHM r = 1.76to). 
Compared to the steady propagation at z < 2.2 cm, 
Fig.[2|a) confirms a deviation of t^{z) from the initial lin- 
ear trajectory by — 600fs at point B (2; = 4.5 cm). This 
temporal shift is accompanied by a 4% change of pulse 
energy cx I2 and by a more than twofold increase of peak 
power Po{z)^ [solid and dashed curves in Fig. [2|b), re- 
spectively]. Pulse duration scales accordingly from an 
initial 20 fs (FWHM) to sub-10 fs at B. Furthermore, 
a Fourier analysis indicates that the center wavelength 
Xo{z) of the soliton shifts from 1060 to 985 nm within the 
2.3 cm propagation from A to B, reflecting the according 
energy transfer. 

Rogue waves, subject to non- Gaussian statistics, have 
been shown to appear in the fiber SC generation both 
with and without Raman frequency shift in [16 . The 
rogue event regarded as an emerging single "champion" 
soliton was linked to multiple collisions between optical 
solitons. As in our case there is no other soliton in reach 
from Ato this parameter change can only be explained 
by nonlinear continuum-soliton interaction. Namely, for 
each soliton velocity there is a spectral slice of disper- 
sive non-solitonic radiation which propagates at nearly 
identical group velocity. Group- velocity matching signif- 
icantly increases the nonlinear interaction length between 
continuum and soliton |2Tl [23] and may lead to the re- 
shaping of the latter [26] . Comparing to the energy of the 
DW that is in temporal overlap with the soliton [dotted 
line in Fig.[2];b)], it is striking that changes of each of the 
soliton parameters t*(z), Ao(^), and Po{z) are strongly 
correlated with the strength of the DW, see positions A 
and A' marked in Fig. [2] It appears still surprising that 
a DW with less than 10% of the soliton amplitude can 
affect its properties so strongly. 

For further investigation of this scenario, we numeri- 
cally isolated the primary soliton and selected segments 
of the DW in Fig.jlja) right at the onset of the trajectory 
curvature, allowing for a deterministic interpretation of 
the acceleration process of the soliton uncoupled from the 
SC generation process. To this end, we inject into the 
fiber a fundamental soliton at As = 1030 nm of 26.6 fs 
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FIG. 3: (a) Time domain evolution along the fiber represent- 
ing the isolated trajectory of a fundamental soliton. Acceler- 
ation results from a cascaded scattering with three DWs. (b) 
Relative change of soliton parameters: (i) photon number (see 
Eq. ^ and 17 ), (ii) energy, (iii) peak power. The deviation 
of a parameter ( is defined as [C(t) — C(0)]/C(0). For clarity, 
deviations (i) and (ii) have been stretched by a factor 10. 
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FIG. 4: (a) Exemplary concave group delay /3i = l^'i^) 
and related group- velocity dispersion /32 = 13" {u)^ with the 
extracted wavelengths for the fundamental soliton at As = 
1030 nm and a dispersive pulse at \d — 614 nm (dashed line). 
Arrows indicate the induced change of {32 and /3i for the soli- 
ton. (b) Spectral evolution along the fiber representing a 
cascaded scattering with three DWs from a soliton. 



FWHM duration together with slightly slower propa- 
gating 53.2 fs time segments of DWs near the velocity- 
matched wavelength of A^^ = 614 nm. 

Figure [3|a) demonstrates an example for a continuum- 
soliton scattering process with three DWs in a suitably 
chosen reference frame. Each collision leads to a stepwise 
acceleration of the soliton, clearly confirming the transfer 
of energy and the concomitant gradual increase in peak 
power [Fig. [sj^b)] as previously seen in Fig. |2] As all ob- 
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scuring continuum components have been eliminated, the 
role of the dispersive radiation can now be seen in much 
greater clarity. The DWs initially propagate at slightly 
lower group velocity than the trailing soliton so that they 
eventually collide. In this collision, the soliton can never 
pass the DW as would be expected in a purely linear- 
optical encounter. Instead, XPM between DW and soli- 
ton causes a frequency shift towards the ZDW (842 nm), 
decreasing the center wavelength of the soliton and shift- 
ing the DW toward longer wavelength. Comparing to the 
underlying dispersion profile [Fig.Qa)], both these shifts 
lead to an acceleration for the respective type of radia- 
tion, as is clearly confirmed by the trajectory curvatures 
in Fig. da). 

Under similar conditions, the impenetrability of the 
soliton trajectory was referred to as an optical event 
horizon for the DW [25l [26]. Its origin lies in a non- 
linearly induced increase of the group velocity caused by 
the leading edge of the soliton. The only way for the DW 
to escape from the event horizon is a shift towards the 
ZDW, i.e., both types of radiation therefore experience 
a strongly enhanced effective XPM. These processes are 
completely elastic, causing a mutual shift of optical fre- 
quencies but never transferring photons from the normal 
dispersion regime into the soliton regime or vice versa. 
With the photon number of the soliton practically con- 
served, the soliton blue shift accordingly causes a mild 
increase of its energy [Fig. [sjb)]. The shifted soliton also 
experiences a considerably smaller ^2 [Fig. |4ja,b)]. Now 
consider that the energy of a fundamental soliton can be 
expressed through Pq and as £^ = 2y^Po|/32|/7. Ob- 
viously, the decrease of cannot be compensated by a 
reduction of E because E also grows. As 7 does not vary 
appreciably, consequently, Pq is forced to grow massively. 
This clearly explains our observations. We repeated these 
simulations with several segments of continuum to prove 
that even higher soliton peak powers can be achieved 
with segments of the continuum containing more energy. 

The nature of the newly observed continuum-soliton 
scattering processes is markedly different from soliton- 
soliton scattering. As continuum radiation quickly dis- 
perses, there will always be temporal slices of the DW 
that effectively copropagate with a given soliton, mak- 
ing mutual extended interaction much more likely than 
the appearance of soliton-soliton processes. We therefore 
suggest that this mechanism contributes to the dramatic 
amplitude increases seen in experimental work [8 . 

The mechanism proposed here does not presuppose 
any special nonlinear effects that are unique to optical 
systems. In comparison to previously discussed mecha- 
nisms of rogue wave formation, our approach essentially 
only presupposes a nonlinear Schrodinger type scenario, 
with a reactive nonlinearity and a concave dispersion pro- 
file, the latter enabling copropagation of radiation with 
opposite signs of dispersion with equal group velocity. 
These conditions are met in a variety of systems, e.g., for 



gravity-capillary waves [27] . Our explanation is therefore 
immediately applicable to a much wider class of physical 
systems. Consequently, we believe that the previously 
disregarded scattering of DWs off solitons opens a new 
perspective on the fascinating appearance of extreme- 
value wave phenomena. 
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